Two non-comoving stiff fluids in radial motion and spherical symmetry 
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The problem of two stiff fluids (energy density — pressure) moving radially in spherical symmetry is 
treated. The metric ansatz is chosen spherically symmetric, conformally static with a multiplicative 
separation of variables. The first fluid is described mathematically via a massless scalar fleld. The 
coordinate system is chosen comoving with the second fluid which the separation of variables requires 
to be stiff too. The fluids are interacting only gravitationally and their energy momentum tensors 
are separately conserved. The Einstein equations are reduced to a single nonlinear ODE of second 
order which is shown to lead to an Abel ODE. A few particular exact solutions were found using a 
polynomial ansatz. The two non-comoving gravitational sources in the solutions can be interpreted 
either as scalar fields or stiff fluids. A complete analysis is performed on the range of parameters 
for which the stiff ffuid interpretation is physically acceptable. General formulas are derived for 
the conformal vectors of the solutions. By making the second fluid vanish, a few single scalar field 
solutions are generated some of which appear to be new. All solutions considered in this paper have 
a time-like singularity at the origin (except the trivial FRW one) and are not asymptotically fiat 
(except the static one with k=0). 

PACS numbers: 04.20.-q 



I. INTRODUCTION 

The most economical way to explain the plethora of 
current cosmological data is the ACDM model which con- 
tains dark matter. Its interactions with the baryonic 
matter is assumed so weak that it still evades a direct 
experimental detection. The various matter components 
in ACDM have the same four-velocity. That is expected 
if they originated from the same component in the past 
or if they have had enough interaction to equalize their 
expansion rates. The existence of a weakly interacting 
component opens the possibility for non-comoving mo- 
tion. Such a scenario cannot be realized within the ho- 
mogeneous and isotropic FRW metrics. Those are di- 
agonal in the usual comoving coordinates and require 
a diagonal energy momentum tensor through the Ein- 
stein equations. Non-comoving components described as 
perfect fluids would necessarily produce off-diagonal el- 
ements in the energy-momentum tensor and are incom- 
patible with the assumed homogeneity and isotropy of 
the ACDM metric. 

Spacetimes with less demanding symmetries, of which 
the simplest are the plane or the spherically symmet- 
ric ones, allow naturally for a non-comoving motion of 
their matter sources. The published exact solutions for 
two non-comoving fluids are scarce due to the prohibitive 
mathematical complexity. The case of plane symmetry 
and two stiff fluids in irrotational motion was solved in 
Shear free two-fluid solutions were considered in 
but the interaction between the fluids was not fixed on 
physical grounds, a problem common to many exact solu- 
tions. Two non-comoving dusts in spherically symmetric 
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radial motion were treated in [J 0] where, by imposing 
the existence of a Killing vector, the problem was re- 
duced to numerical integration of two coupled non-linear 
ODE'S. 

Obtaining non-comoving exact solutions is important 
not only as toy cosmological models but also, from purely 
mathematical point of view, they provide examples of 
solving systems of coupled nonlinear differential equa- 
tions for which general methods yet do not exist. Non- 
comoving fluids could also arise in studies of star evolu- 
tion, for example neutron stars 

When the system of Einstein PDE's is intractable, the 
usual methods to reduce the complexity are separation 
of variables or imposing the existence of Killing or homo- 
thetic vector fields which convert the PDE's into ODE's. 
Stiff fluids (energy density = pressure) appear often in 
the literature since they allow for a simple mathematical 
description in terms of massless scalar fields [gj . 

The present paper was inspired by in which ex- 
act solutions were derived in the context of a modified 
GR theory. As the authors point out, their equations 
are mathematically equivalent in GR to a scalar field 
conveniently interacting with a perfect fluid in spheri- 
cal symmetry. That suggested a similar GR solution is 
also possible for a stiff fluid (described by a scalar field) 
interacting only gravitationally with another fluid. 

Analogously to Q the metric ansatz in the present pa- 
per is chosen spherically symmetric, conformally static 
with a multiplicative separation of variables. The first 
fluid is described by a massless scalar field. The coordi- 
nates are chosen comoving with the second fluid which 
the separation of variables require to be stiff too. The flu- 
ids are interacting only gravitationally and their energy 
momentum tensors are separately conserved. The Ein- 
stein equations are reduced to a single nonlinear ODE of 
second order which is shown to lead ultimately to an Abel 
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ODE. A few particular polynomial solutions are found. 
The two non-comoving matter sources in the solutions 
can be interpreted either as scalar fields or as stiff fluids. 
A complete analysis is performed on the range of param- 
eters for which the stiff fluid interpretation is physically 
acceptable. General formulas are derived for the confor- 
mal vectors of the solutions. By making the second fluid 
vanish, a few single scalar field solutions are generated 
some of which appear to be new. 

The units employed throughout this paper are SttG = 1 
and the metric signature is (- + + +). 



II. CORRESPONDENCE BETWEEN A STIFF 
FLUID AND A MASSLESS SCALAR FIELD 

As was shown in 6] for any perfect fluid with a nonzero 
pressure, the energy momentum tensor and its conserva- 
tion can be written in terms of a scalar function. The stiff 
fluids enjoy particularly simple equations that are math- 
ematically equivalent to those of massless scalar fields as 
outlined below. 

The energy momentum tensor of a stiff fluid in the 
usual form is 



T., = 2eUM, 



(1) 



where e is the energy density (and pressure) and [//^ is the 
four-velocity. If the motion is irrotational, is propor- 
tional to a gradient of a scalar function $, and normalized 



(2) 



where comma denotes partial derivatives. The energy 
momentum conservation, T''^^ = 0, projected parallel 
and perpendicular to leads to the two equations [1]: 



We,, 



(3) 
(4) 



One can always redefine <& ^('i') with an arbitrary 
function F , without changing J7^, until ^ is satisfied 
identically [6]. Then the energy density is given by 



and ([3]) leads to the equation 
1 
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(5) 



(6) 



Substituting ([2]) and ([5]) in (1), obtains the energy mo- 
mentum tensor 



(7) 



which up to a factor of two is the usual tensor for a 
scalar field of zero mass and a vanishing potential. The 
conservation equation of (7) is indeed (6). 

For fluids that are not stiff, the equations analogous to 
(6) and (7) are more complicated and not amenable to 
the separation of variables employed later in this paper. 

III. METRIC AND EINSTEIN EQUATIONS 

The most general diagonalized spherically symmetric 
and conformally static metric with multiplicative separa- 
tion of variables is 



F{t) { -A^{r) + B^{r) dr^ + C\r) dn^ } , 



(8) 

where dfP = dO^ + siv? {6) d(lP' . The above notation is 
similar to the one in [7] but here F(t) is not squared and 
is pulled in front as a conformal factor. That makes the 
equations for F(t) easily solvable. 

The total energy momentum tensor of the two compo- 
nents is 



(9) 



where T^i""* is given by (7) and represents a stiff fluid mov- 
ing along the radial direction; T^fJ represents a perfect 
fluid which will necessarily turn out to be stiff too: 



(10) 



Just like in 7], for simplicity, the coordinate system is 
chosen comoving to the four-velocity 



V = ( 



1 



F^/^\A\ 



,0,0,0) 



(11) 



of the second fluid - the one denoted with "c" for "co- 
moving" . Note that V^^ has zero shear. 

There are four independent non-zero Einstein equa- 
tions G^u = R^u - {Rl2)g^,i, = T^^: 
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(12) 
(13) 
(14) 

(15) 



r 



with dot denoting d/dt and prime denoting d/dr. As 
usual in the separation of variables method, the simpli- 
fying equation is the off-diagonal one (13). Following 
an additive separation of variables is chosen for the scalar 
field 



$(r,i) = W{r)+T(t). 



(16) 



Substituting that in (13) and solving for the scalar field 
functions, W and T, in terms of the metric functions gives 



T(t)-yln 
1 



Fit) 
/Air) 



Wir) = -\n 
Jo V 



(17) 
(18) 



where Tq, Fc, and Ac are constants. The values of Ac and 
Fc are of no importance since they do not influence the 
derivatives of $ which appear in the Einstein equations. 

Note that at the boundary points of the r and t in- 
tervals on which the metric is defined. Air) or Fit) van- 
ish and the scalar field $ becomes infinite. Vanishing of 
Air) signifies a timelike singularity, vanishing of Fit) - a 
spacelike (cosmological) singularity. 



IV. SEPARATE ENERGY MOMENTUM 
CONSERVATION 

The Einstein equations are incomplete without spec- 
ifying the energy transfer between the two fluids. The 
simplest scenario is to assume they interact only grav- 
itationally i.e. their energy momentum tensors are di- 
vergenless. The conservation equation for the comoving 



fluid, T 



[c)fj.i' 



■^1/ — gives the following two equations 
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Assuming a barotropic equation of state, 
ap, a =const, the solution of (19) and (20) is 

Po 



P 



Pir,t) 



F(i)3(i+")/2y4(r)(i+°)/" ' 



(19) 
(20) 

(21) 



where po is an integration constant. The various terms in 
(12) depend either on r or on t. The same should apply 
to the term A^Fp to enable the separation of variables. 
The term is a function of t only if A(r)=const or a = 1 
(stiff fluid) and it is a function of r only if F=const or 
a = —1/3. In the case of A=const, eq. (18) leads to 
$' = which by eq. (2) implies ~ i.e. the two 
fluids are both comoving with the coordinate system - 
that case will be considered further in section IVIIII The 
case of F(t)=const, using eq. (17), leads to $ = which 
by eq. (2) implies the four-velocity U has a zero time 
component which is not physical. The case a = —1/3 
does not describe a known type of matter. The only 
possibility remaining is that the second fluid is also stiff 
with an energy density of 



p(r, t) = p(r, t) 



Po 



F(i)3A(r)2 



(22) 



The conservation equation of the first non-comoving 
stiff fluid, 7"'"''"j'j^ = 0, leads to the equation 



Substituting (17) and (18) in (23) gives 
,A'\' A 



(23) 



C^- 



B 



BC^ 



TlF_Tl 
2F 2 ' 



where 



k = F/F 



(24) 



(25) 



is a separation of variables constant. 

The Einstein equations guarantee the conservation of 
the total energy momentum tensor. Therefore, if one 
of the tensors T''^^^" and T'^"^'"' is conserved, so is the 
other. Thus the Einstein equations and (22) are suffi- 
cient to specify the system; eq. (24) is their consequence. 
Nevertheless, it can and will be used later to simplify the 
analysis. 
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V. EQUIVALENT SYSTEM OF EQUATIONS 

The Einstein equations that need to be satisfied are 
(12), (14), and (15) in which (17), (18) and (22) were 



The first equation above corresponds to 2{GrrA'^/B^ + 
GeeA^/C^), the second to GrrA^/B"^, the third to 
2{Gtt + GrrA^/B^ + 2GeeA^/G'^), all separated in vari- 
ables. The rightmost sides of the equations are the sep- 
aration constants determined in the foUowing way. The 
left hand side of (28) is just twice the conservation equa- 
tion(24). As mentioned before, (24) is redundant and 
should follow automatically from (26) - (28). Neverthe- 
less, it can be used as a shortcut implying that the right 
hand side of (28) is kT^. Another consequence of (24) is 
that F/F = k. That allows to solve (28) for the expres- 
sion: 

fJ(7^o -3) + 4|| = Mro^-3). (29) 

Substituting that in (26) and (27) produces their right- 
most sides. 



VI. SOLUTION OF THE TIME EQUATIONS 

The principal time equation is (25) with the solution 

{/i exp{Vk t) -f /2 exp(-\/A: t), A: > 
ht + f2, k^O (30) 

/isin(^t + /2), fc<0, 

where /i and /2 are constants. In order to satisfy the 
rightmost sides of (26)-(28), F(t) must satisfy (29). Sub- 
stituting (30) in (29) determines the required values of 
Po- 

( k{T§ - 3)/i/2, k > 
Po^l -(r2-3)AV4, k = (31) 
I fc(T2 - 3)AV4, k < 0. 



substituted. A more convenient system of equations, con- 
taining fewer terms, is obtained by forming independent 
linear combinations of those: 



(26) 

(27) 
(28) 



VII. RADIAL GAUGE AND RADIAL 
EQUATIONS 

A coordinate transformation involving only the radial 
coordinate, r = r(f), induces the following transforma- 
tion in the metric functions 

A{r) i(r) = A{r{rj) 

B{r) B{f) = B{r{f))dr/df (32) 

C(V) ^ C{P) = C(r(f)). 

The coordinate freedom in r can be used to select a co- 
ordinate system in which the Einstein equations have a 
simple form. Many choices are possible but a glance at 
(26) and (28) suggests AB = 1 in the new coordinates. 
The same gauge was used in Q- From (32), the required 
transformation, r(f), is obtained by solving the differen- 
tial equation 

A{r)B{r)^^l (33) 
ar 

for f(r) and inverting that to r(f). The integration of 
(33) and the inversion are always possible since A and B, 
being metric functions, should be nonzero. The solution 
f (r) is fixed up to an integration constant implying that 
all radial coordinates for which AB — 1 are related by 
translations: f * — 7~ + const. That residual freedom in f 
will be apparent in the Einstein equations later and can 
be used to absorb constants. From now on, it will be as- 
sumed that the necessary coordinate transformation has 
already been effected and AB = 1 in the new coordinates 
(dropping the tildes for convenience). In this gauge, the 
radial equations (26)- (28) simplify to 

^{{G'Ay-2)=kiTS + l)i3i) 




^VmV 1 o^-^-.|!(T„^-3).4|i.MTo^.l) 

^'(To^-3) + ^ = t(7^o^ + 1) 



I 
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VIII. RADIAL SOLUTION FOR k=0 

When fc = 0, (34) implies 



(37) 



where ri and r2 are integration constants. Substituting 
— R/A^ and fc = in (36) and integrating twice gives 



= Al exp a 



dr 

W) 



(38) 



where a and An are constants. Now C can be found: 



C = exp [ ~a 



dr 



(39) 



The last equation to be satisfied, (35), will put restric- 
tions on a,ri, and r2. Substituting (38) and (39) in (35) 
leads to 



R 



12 



(40) 



Inserting the explicit form of R{r) in the above gives the 
restriction 

r^-r^^ (41) 

An analogous solution was obtained in Q ■ The constant 
a must be real since it enters the derivative A' /A ~ 
a/{2R). This is possible only if ri and r2 in (41) have 
equal imaginary parts. The roots ri and r2 must be com- 
plex conjugate to each other to produce the real-valued 
polynomial R{r). The last two statements are compati- 
ble only when ri and r2 are real. The two distinct cases, 
ri = r2 and ri ^ r2 are now discussed. 

The case A = const of two comoving fluids was men- 
tioned in section ITVl It implies fc = from (36), a — 
from (38) and ri = r2 from (41). The metric is spatially 
flat FRW with Big Bang or Big Crunch depending on the 
sign of To as will be discussed later. The energy density 
(5) is non-negative if A = const. The energy density (22) 
of the other fluid is non-negative if Tq < 3, see (31). 

In the case of ri ^ r2 the two fluids are not comoving 
with each other. The metric functions become 



A\r) 



- 42 

- ^0 



r2 



(42) 



-^(r-ri)(r- 



r2) 



r2 



where the power n = ±(1 + 1 / Tq) The sign of n is 
the sign in (41) chosen for a. 



IX. PARTICULAR RADIAL SOLUTIONS FOR 

/ 

For nonzero k, the radial equations (34)-(36) can be 
reduced to a single equation involving the function 



S{r) = A^C^, S" + 2. 



(43) 



The heuristics behind taking that combination is trying 
to imitate the way the fc = solution was obtained as 
far as possible. A big advantage of that choice is that 
equations (34) and (43) can be solved algebraically for 
^2 and 6*2 in terms of 5*: 



A' 



C = 







S" 


- 2 


S{S" 


-2) 


k{Ti- 


f 1)' 



(44) 



where S > because of (43) and k/{S" — 2) > because 
of (34). Substituting the above in (36) obtains 



S' 



SS" 



2T2 



(5" -2). 



(45) 



S"-2J + 1 

Integrating that equation once gives 

(ro2 + i)^5"' ^ {S"~2){S'{l-T^) + 'irTS + Ci), (46) 

where Ci is the integration constant. Equation (46) can 
be integrated again resulting in 

(T2 + 1)SS" ~ {S')^ - S"(4T2r + C\) + 2{T^ + 1)5* (47) 

+ (4r2r2 + 2Cir + C2) = 0, 

where C2 is the new integration constant. The equation 
left to be satisfled is (35) - it will put restrictions on (47). 
Substituting (44) in (35) gives 



S"'S 
S"'-2 



+ 2(ro2 + i)s" 



-(3To2 - 1)S"2 -4r(2(S"' 



S"'S 
S" -2 
2)S = {) 



(48) 



The order of this equation can be reduced by solving (46) 
for S"'S/[S"-2) and substituting that back in (48). The 
result is 

[T^ -f \)SS" - - S'{AT^r + d) + 2(To2 + 1)5* (49) 

^ " Cif = 0. 



Comparing that to (47) reveals that the restriction put 
by (35) on (47) is simply C2 = Cl/{AT^). 

In the above treatment, the Einstein equations (34)- 
(36) were reduced to (49) - a single non-linear ODE of 
second order. The integration constant Ci in (49) ex- 
presses the translational freedom of the radial coordinate 
discussed in section IVlIl if S(r) is a solution of (49) for 
a given constant Ci, then S{r) — S{r — ro),ro — const 
is also a solution of (49) but for another constant Ci = 
Ci - iT^ro. 

Three particular solutions of (49) were found by trying 
a polynomial ansatz for S(r). These are: 



Sir) = 
A'ir) - 



2Tr 



rp2 

-^0 



"-^{r-ro)' 



(50) 



Vk \To{7 



To 



T2 ■ 
^0 



^o)| 
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where Tq > 1, r 7^ tq, k > 0; 



S{t) = d^ir^r^f 



(51) 



where Tq 



A^{r) = ^i-k/3){d-{r-ror) 
C\r) - ^{-i/k){d-{r-r^Y), 

= 1/3, d > 0, fc < 0, |r - rol < Vd; 



S'(r) = a(r - ro)'' + (r - ro)^ + fe(r - ro) 



(52) 



3 \/ a V r - To 



C2(r) 



where 1/3, a 7^ 0, S'(r) > 0, fc/a > 0. 

In the above solutions, a, 6, d are constants and tq is 
a constant reflecting the translational freedom in r. The 
constraints after each solution come from the require- 
ments S* > and k/{S" — 2) > mentioned after equa- 
tion (44) and from requiring positivity of the expressions 
inside the square roots . Since the solution (52) depends 
on two constants, a and 5, it captures the a ^ part of 
a general solution for Tq = 1/3. 

Although the main equation (49) is not particularly 
complicated in appearance, its general solution for arbi- 
trary To could not be found. The following sequence of 
variable transformations: Ci = (by translation in r); 
S{r) = u{r)r'^; r = exp(a;); reduce (49) to an equation 
for u{x) in which the independent variable x does not 
appear explicitly: 



(T^ -I- l)uu^., + (3r2 - l)(u, -2)u-z 
+2(T2 - l)w2 



0, 



where subscripts denote ordinary derivatives. The stan- 
dard technique for solving such type of equations is the 
substitution for the first derivative = q{u) (hence 
Uxx — <l{u)qu) which converts the above equation to an 
Abel type ODE for q{u): 



qqu 



{-u - ATg + STo^u) q 



u{TS + l] 
2(u- 



uiT§ + 1) 



Tiu- 



u{T§ + l) 



0, 



for which general methods of solving are not known at 
present. Two particular solutions of (49) are known, 
which translate into particular solutions of the above 
Abel equation. These are solution (50) and the following 

S{r) = [r - rof + b{r - ro) + 



4(T2 



1 



The above S{r) was not mentioned before since it has the 
forbidden property S" = 2 and thus it cannot generate 



solutions for and C^. For Riccati types of ODE's, 
which are particular cases of Abel, there are techniques 
to obtain general solutions starting from one or more 
particular solutions ^ . Unfortunately, the equation here 
is of the more general Abel type for which such methods 
are not currently available. 



X. PHYSICALITY OF THE STIFF FLUID 
INTERPRETATION 

In any solution, the metric functions FjA^, and 
are required to be positive to produce the right metric 
signature. Additionally, if the two matter sources are in- 
terpreted as stiff fluids, their energy densities must be 
non-negative and the four-velocities must be future ori- 
ented. The energy density of the comoving fluid is non- 
negative when po > 0. The energy density (5) of the first 
fluid written more explicitly is 



4FA2 



(A2 



To 



(53) 



and its four-velocity (2) is 



TqF (A^ 



2FVe V A'^F' To 



0, 



(54) 



The four- velocity is future oriented, J7* > 0, when TqF < 
0. Hence the coordinate system must be expanding, F > 
0, when Tq < and contracting, F < 0, when Tq > 0. 
The allowed range of the time coordinate is restricted by 
the inequalities F>0,F<Oor >0 (depending on the 
sign of Tq) and the non-negativity of e: 



> 



Tq' 



(55) 



The allowed intervals for r are determined by the con- 
ditions A^ > 0, > and are bounded by points at 
which A^ or become zero or infinite. 

The fc = solutions in (42) have a time function F{t) = 
fit + f2- If To < (expansion), the allowed time interval 
is t > — /2//1, /i > to satisfy the conditions F > 
0, T > 0. At time t — /2//1 the energy densities (22) 
and (53) become infinite signifying a spacelike singularity, 
a Big Bang. Note that the velocity component U"^ can be 
negative if the first fluid is expanding slower than the one 
comoving with the coordinates. Analogously, the To > 
case (contraction) describes a spacetime ending with a 
Big Crunch at t ^ — /2//1 and the time is restricted to 
t < — /2//1, /i < 0. If r2 denotes the larger root in (42), 
positivity of the metric function requires r < ri or 
r > r2- Without loss of generality, one can take for the 
allowed range r > r2. A simple translation in r can set 
r2 to the conventional zero value for the origin. Since 
\n\ < 1, C^ir) does not become infinite even at ri or 
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r2 . If n > (n < 0) the function (r) vanishes at 
r — r2 {r — ri) which makes the energy densities (22) and 
(53) infinite, signifying a timeUke singularity at r2(ri). 
The energy density (22) is non-negative if Tq < 3, see 
(31). The hardest physical requirement to satisfy is the 
non-negativity of e, (55). Written explicitly, this is 



fi 



^ Vg(r^ + i) 



5F2 - n\ 



ri|i+"|r- 



r2 



(56) 



Since \n\ < 1 and ri ^ r2 in the non-comoving case, the 
above inequality restricts r for fixed t or restricts t for 
fixed r. Thus the energy density e is not non-negative 
over the full allowed ranges of the coordinates. A similar 
problem occurs sometimes for a single stiff fluid in plane 
symmetry |f3]. There the energy momentum tensor was 
reinterpreted as an anisotropic fluid in regions of negative 
e. If the non-comoving component in the present article 
is interpreted simply as a scalar field instead of a stiff 
fluid, the problem of negative e becomes irrelevant. 

The fc > solution (50) has a time function F ~ 
fi exp{-</k t) -t- /2 exp(— Vfc t). When the coordinates are 
expanding (Tq < — 1 < 0), the conditions F > (metric 
function) and F > Q (future oriented four-velocity) lead 
to the restrictions 



/i > 0, i X, 



1 



2^/k 



In 



(57) 



The corresponding restrictions for contracting coordi- 
nates (To > 1 > 0) arc 



/2 > 0, t<t,= 



1 



2^/k 



In 



(58) 



The non-negativity of e condition (55), written explicitly 
is 



/i exp(%/fc t) - /2 exp(-Vfc t) 
fi exp(\/fc t) -I- /2 exp(-Vfc t) 



> 



iTn 



(59) 



For expanding coordinates (Tq < — 1, /i > 0), (59) 
translates into 



exp 



{2Vk t) > 



|To| + 1 h 
1 ■ h' 



Tn 



(60) 



If /2 < 0, (60) is not restrictive at all and the final allowed 
range of t is given by (57); the spacetime starts with a 
Big Bang at t*. If /2 > 0, (60) is more restrictive than 
(57) and the allowed range of t is given by 



t > 



1 



2Vk 



In 



ITn 



Tn -1 



h 

h 



(61) 



In this case the stiff fluid interpretation fails at times 
too close to the Big Bang time . The case of contract- 
ing coordinates (To > 1, /2 > 0) has a similar analysis 
resulting in the following allowed ranges: 



t < i,, /i < 
1 . ^To 



(62) 



t < 



2^/k 



In 



Tn 



/i >0. 



When /i > the stiff fluid interpretation fails at times 
too close to the Big Crunch time The metric func- 
tions in solution (50) vanish at r = ro. Without loss of 
generality, the allowed range can be taken as r > ro. A 
translation of the r coordinate can always set vq to the 
conventional zero value for the origin. 

The A: < solution (51) does not allow for a stiff fluid 
interpretation over the full range of t set by the T > 
condition. Similarly to the k = solution (42) discussed 
before, the energy density e becomes negative when r 
approaches the roots tq ± Vd. The problem stems from 
{A^y becoming infinite at the roots. The detailed anal- 
ysis is left to the reader. 

Solution (52) allows for a stiff fluid interpretation when 
the derivative 



2\/ 



1 [k a(r-ro)-6/(2(r-ro)^) 
3 V a ^a(r-ro)2 + l + 6/(r-ro) 



(63) 



remains finite at the end points of the allowed r intervals. 
The end points are roots of S{r) which are also zeros of 
C^(r). The denominator in (63) vanishes at any root of 
S{r) except tq. The only way for (63) to remain finite at 
such a root, r» ro, is that its numerator also vanishes 
leading to the system of equations 



2(n - ro)2 







(64) 



air* — ro) + 1 + 



= 0. 



r* - ro 



The solutions of the system are b — ±2/(3V— 3a), = 
ro=F 1/v— 3a, where a < for real values. Unfortunately, 
S{r) < in the neighborhood of r* hence r^, cannot be 
an edge point of an allowed r interval and is of no fur- 
ther interest. The root tq has a different behavior. The 
derivative (63) is finite at r = ro only when 6 = 0. If 
b — and a < then (63) is finite at ro but infinite 
at the other two roots ro ± Xj \J—a and the stiff fiuid 
interpretation fails close enough to the roots. 

The case 5 = and a > (A: > 0) is more benign since 
then S'(r) is positive everywhere except at the double root 
r = ro. The allowed r interval can be taken as r > ro. 
The derivative (63) is zero at ro and increases asymp- 
totically towards the value VA;/3. The non- negativity 
condition (55) will be satisfied over the whole allowed r 
interval if |T/T| > 'i\{A^)'U^^ = ^fk: 



fi exp(\/fc t) — /2 exp(-\/fc t) 
fi exp(%/fc t) + /2 exp(-Vfc t) 



> 1. 



(65) 



When To = — 1/ \/3 < (expansion) , the conditions F > 
0, F > again lead to the restrictions (57). Taking into 
account T > 0, T > to open the absolute value in (65) 
leads to /2 < so the final result is 



To 



4/73, /i > 0, /2 < 0, i > U 



(66) 
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This spacetime has expanding coordinates and a comov- 
ing fluid that starts with a Big Bang at . 

In the opposite case To = > (contraction), 

the conditions F > 0, F < again lead to the restrictions 
(58). Using F > 0, F < to open the absolute value in 
(65) leads to /i < and the final restriction on t is 



To - +1/V3, /i < 0, /2 > 0, t < U 



(67) 



This describes a spacetime with contracting coordinates 
and a comoving fluid that ends in a Big Crunch at . 



XI. CONFORMAL VECTOR FIELDS 

It is important to know if the spacetimes obtained 
posses Killing, homothetic or more generally conformal 
vectors. By definition, a conformal vector field, sat- 
isfies 



-2? 9f^i^ = -^(x) gi,v 



(68) 



where £ is the Lie derivative along ^ and "0 is a scalar 
function. Here g^i, is the metric (8) in the B = 1/A 
gauge. For simplicity, the conformal vector is restricted 
to the form ^'^ — (a, 6, m, n) where its components are 
functions of (i, r) only. The off-diagonal components of 
(68) give m — const, n = const and 

b = A^a'. (69) 
The diagonal components of (68) lead to to = and 

'A' C 



b' = b 
a = b 



A G 
C^_A^ 



^{r,t) = a{t,r)j+2b{t,r)^. 



(70) 
(71) 
(72) 



Equation (70) can be integrated leading to the equivalent 
set of equations 



b = f{t)AC 

a = f{t){AC' ~CA') 



(73) 
(74) 

(75) 



0(r,t) = a{t,r)j+2b{t,r)^, (76) 

where f(t) is an arbitrary function of time. 

The solutions of (73)- (76) split into two classes. In the 
case of f{t) = 0, the solution is 

= (a, 0, 0, n), ^ = aF/F, (77) 
where a and n are arbitrary constants. 



Additional solutions for f{t) ^ are possible only if 
the consistency condition (a)' = (a')' : 



^{AC" -CA") ^ const 



(78) 



is satisfied. That condition (with const=0) is true for 
solution (42) with ri = r2 (the FRW case). For solutions 
(50) and (51), the combination AC" — CA" vanishes too 
since C (x A. In all those cases AC" — CA" = and then 
the extra solution for is 

ait,r) = {AC'-CA') J fit)dt + f'{t) J ^dr(79) 

b{t,r) = fit) AC 
m{t,r) = 
n{t,r) = const 

ij{r,t) = a{t,r)^+2b{t,r)^ 

where f{t) is an arbitrary linear function, f'{t) — const. 
Note that the above conformal vector applies only when 
AC - C^'=const. 



XII. SINGLE SCALAR FIELD SOLUTIONS 

By a proper choice of the parameters, the energy den- 
sity (31) of the comoving fluid can be set to zero leaving 
a solution with a single non-comoving stiff fluid which 
can be interpreted as a massless scalar field. The case of 
fc = 0, Tq = 3, n = ±^/3/2 is the non-static scalar field 
solution described in [ly|. The case of /c = 0, /i = 0, Tq 
arbitrary and n arbitrary gives a particular class of 
asymptotically flat static solutions published in [ll[. 

Other possibilities for generating a single scalar field 
solution are fc > 0, /1/2 = 0; fc > 0, = 3; fc < 
0, Tq = 3. These do not appear to have been published 
before. 



XIII. CONCLUSIONS 

A few particular exact solutions were derived for two 
stiff fiuids (or massless scalar fields) in non-comoving ra- 
dial motion in spherical symmetry. The simplest possible 
physical assumption is investigated when the energy mo- 
mentum tensors of the fluids are separately conserved. 
Although that seems too idealized, the same assumption 
is applied to the different components in ACDM. Spher- 
ical symmetry allows for a more complicated thermody- 
namics though. There could be a temperature gradient 
along the distinguished radial direction and correspond- 
ingly a heat flow. The current paper assumes that the 
heat conduction of the fluids is negligible hence they can 
be described as perfect fiuids with no heat flow. Two 
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of the presented solutions can be interpreted as stiff flu- 
ids with a non-negative energy density over the whole al- 
lowed range of coordinates. The other solutions can prob- 
ably be matched as parts of bigger piece-wise solutions. 
Due to the separation of variables anzats (16), all solu- 
tions considered in this paper have a timelike singularity 
at the origin of r (except the trivial FRW one) where A{r) 
vanishes and and a spacelike (cosmological) singularity at 
the origin of time where F{t) vanishes. Moreover, the so- 
lutions are not asymptotically flat (except the static one 



with k=0). 
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